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^s^j ■ Abstract 

We calculate the survival probability of J/Sfr particles in various col- 
liding systems using a Glauber model. An analysis of recent dataffl] has 
reported a J/^-nucleon breakup cross section of 6.2±0.7 mb|| derived 
from an exponential fit to the ratio of J/$ to Drell-Yan yields as a func- 
tion of a simple, linearly-averaged mean path length (L) through the nu- 
clear medium. Our calculations indicate that, due to the nature of the 
calculation, this approach yields an apparent breakup cross section which 
is systematically lower than the actual value. 



1 Introduction 



X 

■ Recently, data from the NA50 collaboration on J/ "J yields from p+ A and A+ A 

collisions have become available, which have been interpreted as displaying 
"anomalous J/^ suppression" H . It has been predicted that in the case of 
a phase transition to a quark-gluon plasma (QGP), the yield of J/ty is sup- 
pressed due to Debye screening S. There have been numerous efforts to explain 
this data, invoking such effects as interactions with co- moving hadronsQ and 
suppressed production due to initial state energy loss(|, |(|, as well as decreased 
yields in a QGP §. 

However, even in the absence of a deconfincmcnt transition, the J/^f can 
undergo inelastic interactions as it traverses the nuclear medium (as has been 
shown recently in data from p+A collisions ||), which reduces the overall J/^ 
yield. Thus, the signature for a deconfinement transition would be J/^S sup- 
pression at a level beyond that expected from the normal nuclear medium, an 
effect referred to as "anomalous J j^l suppression" [p] . 
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In the analysis of the NA50 data, the ratio of observed J/*f> to Drell-Yan 
yields is interpreted as proportional to the J/\E r survival probability. This ra- 
tio is then fit as an exponential function of a simple, linearly-averaged mean 
path length through the nuclear medium (L) which returns a value for the 
J/^-nucleon breakup cross section of 6.2±0.7 mb|^]. The comparison of this 
experimentally determined cross section and theoretical modeling is important 
in determining the relative contributions of the "pre-resonance" color octet and 
the color singlet states ||, |ll|, [l2| to the overall absorption in A- A collisions, 
as well as in estimating the magnitude of the anomalous J/^> suppression. 

In this paper we describe a calculation of J r /\E f survival probability in A+A 
collisions. The results of this calculation indicate that a breakup cross section 
derived in the manner of |l| is systematically lower than the actual value. We 
will discuss the origin of this discrepancy. 



2 The Geometric Model 

In order to study the effect of collision geometry on the J/^P survival proba- 
bility, we have utilized a Glauber model approach, implementing Woods-Saxon 
paramctcrizations of nucleon distributions within the colliding nuclei. Thus, 
assuming spherical symmetry for the nuclei, the nucleon density is described 
by: 



dN 
— oc 
dR 



1 + exp 



R-C 
D 



(1) 



where 



C = Ro x Ay u 3 clcus (2) 

where Ro= 1.11 fm, A nuc i eus is the atomic number of the nucleus, and D = 0.75 
fm. For the specific case of a Pb nucleus, the parameters used are Ro = 1.1 
fm, C = 6.624 fm, and D = 0.549 fm, as given in jl3|. The nucleon-nucleon 
cross section used in the Glauber calculation is Ci„t=30 mb. The results of the 
calculation are not sensitive to the exact value of this parameter over the range 
25-40 mb. 

Production of J/^f is modeled by assuming that each nucleon-nucleon colli- 
sion has the same probability of producing a J/\&. We note that this modeling 
assumes that initial state energy loss does not affect J/\P production. We are 
currently investigating the validity of this condition; the results of this study 
will be addressed in detail in a follow-up paper |3J. 

After production, the J/\P is taken to be at rest in the center-of-mass frame, 
such that the survival probability is a function of the path length through nuclear 
material which the J must traverse, the nucleon density and the breakup 
cross section (J,j/\^-n- Specifically, the survival probability calculation requires 
an integration over z: 

Pltr 3 = \ eM-p(R) vj/v-nz) dz (3) 
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where z tg t,proj is the position of the trailing edge of the target and projectile 
nuclei, respectively (defined by the z value, at the same x and y values as the 
production point, where the nucleon density is negligibly small) and p(R) is the 
local nucleon density. The total survival probability is just the product of the 
two terms, P surv = P*^ x Pj%%. 

Clearly, the survival probability of a given J/^Sf depends strongly on where 
in the interaction region it is produced. However, for any given collision ge- 
ometry (i.e. impact parameter), there are a wide range of possible production 
points, each with a different J/^l survival probability. Thus, the mean survival 
probability for a given geometry represents an average over all possible J/\I/ 
production points within the interaction region. We note that in the inter- 
pretation of actual experimental data, the situation is further complicated by 
the impossibility of experimentally selecting on impact parameter, given finite 
detector resolution and the natural dispersion in the various centrality-related 
measurable quantities. Thus, measured J/vI/ yields represent an average over 
various collision geometries in addition to the average over survival probabilities 
for a fixed, "ideal" geometry, as discussed here. 

In previous analyses]]!], the J/\P survival probability was taken as propor- 
tional to the ratio of observed J r /\E f yields to Drell-Yan. For various classes of 
event geometries, this ratio was fit as an exponential function of (L), the mean 
path length through the nuclear medium which the J/\& must traverse. In that 
approach, (L) for a given class of geometries was the average of L(b) for the 
impact parameters b contained in that class and 

L{b) = j Q (j pd z y (4) 

where po is the mean nucleon density, p is the local nucleon density and the 
integration is performed over the interaction region. We have performed a sim- 
ilar calculation within our Glauber model. In our model, we have used a mean 
nuclear density of po=0-1793 nucleons/fm 3 ; the results of the calculation are 
insensitive to this normalization factor. For a given colliding system, the value 
(L) averaged over all impact parameters (i.e. minimum bias) can be calcu- 
lated analytically (see for details), which provides a cross check of the (L) 
calculation; our values agree well with the expected values. 

Shown in Fig. |l| arc the calculated J/vP survival probabilities for various 
colliding systems and geometries, plotted as a function of (L). To arrive at 
the specific (L) points given in the figure, we have simulated NA50 transverse 
energy binsH; transverse energy is assumed to scale as the number of wounded 
nucleons[[l4|. However, the substance of the figure does not depend on the 
specific centrality bins chosen. For the points shown, a breakup cross section of 
&J/&-N = 6.2 mb is used as input to our Glauber calculation. As can be seen 
clearly in Fig. [3, the apparent breakup cross section of 5.5 mb obtained from 
an exponential fit versus (L) differs substantially from the actual breakup cross 
section. Further, the points from the calculation, while they give a reasonable 
fit, do not lie precisely on an exponential. Since the points do not lie on a perfect 
exponential, the amount of error incurred by fitting as a function of (L) will 
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depend on the exact centrality binning implemented. However, we emphasize 
that the discrepancy between the apparent and actual breakup cross sections is 
inherent in the collision geometry, and would be present even if it were possible 
experimentally to precisely determine the event impact parameter. 

It is important to gauge the level of the deviation from exponential simply 
due to the error in calculating (L), since the anomalous suppression in Pb+Pb 
collisions is measured relative to an extrapolation of the exponential fit. We have 
found that there is a systematic deviation from exponential, but of a relatively 
small magnitude. And it is notable that this deviation for the most central Pb 
+ Pb values is above the projected exponential fit, in the opposite direction of 
the suppression. Thus the observed suppression is not explained away by the 
error in plotting versus (L). 

However, the breakup cross section calculated this way is incorrect, and this 
error is straightforward to understand. For any single J/vI 7 , the survival proba- 
bility follows an exponential in L for its given production point. However, when 
calculating the mean survival probability, care must be taken on how the mean 
path length is calculated. While the calculation of (L) as described above treats 
all points in the interaction region as contributing equally to production, it does 
not take into account the fact that, due to the exponential nature of absorption, 
all points in the interaction region do not contribute equally to net J yields. 
In order to accurately determine the actual breakup cross section from measured 
J /if? yields, the possible production points must be weighted by their survival 
probability, which depends on the breakup cross section; thus, a proper calcu- 
lation of (Jj/^-n could be accomplished via an iterative calculation, similar in 
spirit to the one detailed here, which exactly returns the input cross section if a 
weighted mean path length is used. Since surviving J/vff particles preferentially 
come from near the outer edges of the interaction region, an exponential fit 
versus a linearly-averaged (L) systematically returns an apparent breakup cross 
section which is less than the actual value. In Fig. ^, apparent breakup cross 
sections for various sample input breakup cross sections are shown; a fit to the 
points yields: 

Tactual = 0.118 + 0.932 Apparent + 0.030 al ppaIcnt (5) 

The exact values for this fit correspond to the centrality bins chosen, and vary 
slightly for different binning. With this caveat in mind, for an apparent breakup 
cross section of 6.2±0.7 mb, as reported in 0, our model indicates that the 
actual breakup cross section is 7.1±0.9 mb. 



3 Conclusions 

We have calculated the survival probability of J/ty particles in various colliding 
systems using a Glauber model, and compared the results to the nuclear ab- 
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sorption observed by the NA50 collaboration [Q. While our calculation does not 
reconcile the issue of whether anomalous J/*f> suppression has been observed in 
the NA50 Pb+Pb data, it does indicate that the apparent breakup cross sec- 
tions determined by a fit of yields versus (L) are systematically lower than the 
actual value, and a proper evaluation of the breakup cross section requires an 
iterative calculation. 
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Figure 1: Simulated J/^ survival probability as a function of (L), the mean 
path length through the nuclear medium, for various colliding systems (p+Pb 
are open squares, S+U are solid triangles and Pb4-Pb are solid circles-see the 
text for details on the values of (L)). The points were generated using a breakup 
cross section of 6.2 mb in the nuclear medium. A fit to the calculated survival 
probabilities, shown as a dotted line, returns an effective absorption cross section 
of 5.5 mb. For comparison, the dashed line indicates the trend the points should 
take if (L) were a valid indicator of the mean free path in the nuclear medium. 
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Figure 2: Breakup cross sections as determined by an exponential fit of J/\l/ 
yields as a function of {L), compared to the actual breakup cross sections used 
in the calculation. The value returned from the fit, shown as solid circles for 
several input values of breakup cross section, are consistently lower that the 
actual input value (for comparison, a solid line of unit slope is also shown). 
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